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Abstract 

A field-theoretic description of the critical behavior of weakly disordered sys- 
tems with a p-component order parameter is given. For systems of an arbitrary 
dimension in the range from three to four, a renormalization group analysis of the 
effective replica Hamiltonian of the model with an interaction potential without 
replica symmetry is given in the two-loop approximation. For the case of the one- 
step replica symmetry breaking, fixed points of the renormalization group equations 
are found using the Pade-Borel summing technique. For every value p, the threshold 
dimensions of the system that separate the regions of different types of the critical 
behavior are found by analyzing those fixed points. Specific features of the critical 
behavior determined by the replica symmetry breaking are described. The results 
are compared with those obtained by the e-expansion and the scope of the method 
applicability is determined. 

1 Introduction 

When the renormalization group approach is applied to describe the critical behavior 
of disordered systems with quenched disorder, the method of replicas [H E] is used to 
restore the translation symmetry of the effective Hamiltonian describing the interaction of 
fluctuations. However, in the studies jHEHEj, it was conjectured that the replica symmetry 
could be broken in systems with quenched disorder. In jHEJE], the physical concept of the 
occurrence of numerous local energy minima in disordered systems with random transition 
temperature was used to give a renormalization group description of the 4 model with 
the interaction potential characterized by the broken replica symmetry. For this purpose, 
the e-expansion technique was used in the lower order of the theory. For systems in which 
the number of components of the order parameter p is less than four, it was discovered 
that the breaking of replica symmetry is the crucial factor in the critical behavior. It was 
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shown that, for p in the range from one through four, two modes of the system behavior 
are possible. The first one determines a nonuniversal critical behavior, which depends on 
seed values of the model parameters and, ultimately, on the concentration of impurities in 
the system. The second mode is characterized by the absence of a stable critical behavior, 
as also is the most interesting case of Ising (p = 1) systems. Even though the implications 
of these studies are very interesting, the results of a field-theoretic description of certain 
homogeneous and disordered systems in the two-loop and higher order approximations 
based on the asymptotic series summation techniques [Zj showed that the stability analysis 
of various types of critical behavior that uses the first-order terms of the e-expansion can 
be considered only as a coarse estimate, especially for multivertex statistical models [Hj. 
For this reason, the results of investigation of the replica symmetry breaking (RSB) effects 
obtained in |3J |3 require revaluation from the viewpoint of a more accurate approach. 

To this end, we proposed in [5J in the framework of the field-theoretic approach, 
a renormalization group description of the model of weakly disordered three- and two- 
dimensional systems with the fourth-order interaction potential with respect to the order 
parameter fluctuations, which determines the replica symmetry breaking. An analysis 
of solutions to the renormalization group equations carried out in the two-loop approxi- 
mation with the sequential application of the summation technique for Pade-Borel series 
showed that the critical behavior of three- and two-dimensional systems is stable with re- 
spect to the relative influence of the RSB effects, and the former scenario of the quenched 
disorder influence on the critical behavior is realized [TT] . 

However, the scope of the results obtained in |3J E] remains unclear. In particular, it 
is interesting to establish the threshold dimensions of the disordered system, d c (p), that 
separate the domain of influence of RSB effects from the critical behavior domains in which 
these effects are insignificant. It is also interesting to apply the renormalization group 
approach to analyze the behavior of systems with replica symmetry breaking in which 
no stable critical behavior exists and the strong coupling mode occurs (see |H E]). A 
theoretical analysis of this phenomenon is especially important from the viewpoint of the 
possible manifestation of RSB effects in strongly disordered systems and their observation 
in computer models of the critical behavior under the impurity concentration exceeding 
the impurity percolation threshold when extended impurity structures are formed in the 
system [T2~| . 
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This paper is devoted to the consideration of the above-mentioned problems. For 
weakly disordered systems of an arbitrary dimension in the range from three to four, an 
analysis of the critical behavior of the model with an RSB potential is carried out based 
on the renormalization group approach in the two-loop approximation with the use of 
summation techniques. Our analysis does not rely on the e-expansion technique. 

2 Definition of the model and the calculation proce- 
dure 

The model Ginsburg-Landau Hamiltonian, which describes the behavior of a ^component 
spin system with weakly quenched disorder near the critical point has the form 




where the random phase transition temperature has the Gaussian distribution 5t(x) with 
the variance << (5t(x)) 2 >>~ u, which is determined by a positive constant u and 
is proportional to the concentration of the structure defects. The application of the 
conventional replica method (see, for example, makes it possible to average over the 
temperature fluctuations Sr(x) and reduce the problem of the statistical description of a 
weakly disordered system to the problem of the statistical description of a homogeneous 
system with the effective Hamiltonian 

„ i P n i p n i p n 

H n = J ^2EE[v0-(z)] 2 + 2 r EE[^)F + 4 E gaA€^)} 2 [<P b M)} 2 }, (2) 

i=l a=l i=l a=l i,j=l a, 6=1 

Here, the index a enumerates replicas (images) of the original homogeneous component 
in Hamiltonian and the additional vertex u, which occurs in the interaction matrix 
gab = g^ab — u, , specifies the effective interaction of fluctuations of the (n xp) -component 
order parameter through ground state of the system with the configuration <p{x) = ( 
at T > T c ), is performed at the scale of the correlation length, which turns to infinity 
the defect field. This statistical model is thermodynamically equivalent to the original 
disordered model in the limit n — > 0. 

The subsequent renormalization group procedure, which statistically takes into ac- 
count the contribution of long-wave fluctuations of the order parameter relative to the 
at the transition temperature T c . This procedure makes it possible to analyze possible 
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types of the critical behavior and conditions of their realization and calculate the critical 
indexes. 

However, it was shown in (31 El E] that a macroscopically large number of spatial 
regions with <p{x) 7^ appears in the system due to fluctuations of the random transition 
temperature at [r — 8t(x)\ < 0. These regions are separated from the ground state by 
potential barriers. To describe the statistical properties of systems with multiple local 
energy minima, the replica symmetry breaking formalism (suggested by Parisi) was used 
in pfl El E] by analogy with spin glasses |SJ. According to the reasoning presented in 
jU El E], the statistical calculation of the contribution of nonperturbation degrees of 
freedom associated with the order parameter fluctuations relative to the configurations of 
the field 4>(x) at the local energy minima results (when the replica procedure is applied for 
the weak disorder) in the appearance of additional interactions of the form in Y^ a ,b 9ab ( t ) a ( t ) 'b 
the effective replica Hamiltonian. Here, the final matrix is no longer replica symmetric 
with g a t = g5 a b — u, but rather has the RSB Parisi replica structure [13]. According to 
|H El El E], the matrix g a t with the RSB structure is parameterized (in the limit n — > 
) in terms of its diagonal elements g and the off-diagonal function g(x), which is defined 
on the interval < x < 1: g a b — » (g,g(x))- Here, operations with the matrices g a b are 
defined by the rules 



9ab (9 k ; 9 k (x)), (g 2 )ab = 9ac9cb ~> (5j c(x)), (f) ab = 9ac9cdgdb ~> d(x)), (3) 

c=l c,d=l 

where 

c = f-\ dxg 2 (x), c(x)=2[g- dyg(y)}g(x) - dy[g(x) - g(y)] 2 } (4) 
Jo Jo Jo 

d = eg- dxc(x)g(x), d(x) = [g - dyg(y)]c(x) + [c - dyc(y)]g(x) - 
Jo Jo Jo 

- / dy[g{x) -g(y))[c(x) - c(y)). 
Jo 

The replica symmetric situation corresponds to g(x) = const, which is independent of x. 

The renormalization group description of the model specified by the replica Hamilto- 
nian (j2J was carried out in the framework of the field-theoretic approach in the two-loop 
approximation for systems of an arbitrary dimension in the range from three to four. 
Possible types of critical behavior and their stability in the fluctuation domain are deter- 
mined by the renormalization group equation for the coefficients of the matrix g^. They 
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were determined by the conventional method based on the Feynman diagram technique 
for the vertex parts of the irreducible Grin functions and the renormalization procedure. 
For example, in the two-loop approximation, the two-point vertex function T^, the four- 
point vertex functions T^, and the two-point function T^' 1 ^ with the insertion (0") 2 have 



the form 






dk 2 




r (4) 
1 ab 





= l + Afg 2 aa + 2pfJ29ac9ca, (5) 

c=l 

n 

= 9ab~pJ2 9ac9cb ~ ^9aa9ab - Agfa, + (8 + 16%^, + (24 + 8h)g 2 aa g ab + 

c=l 

n n n 

+ A8hg aa g 2 ab + Ag aa g bb g ab + 8ph ^ g ac g 2 cb + 8phg ab ^ g ac g cb + Aphg ab ^ g 2 ac + 

c=l c=l c=l 

n n n 

+ 2p g ac 9cc9cb + 4pg aa J2 9ac9cb + p 2 J2 9ac9cd9db, (6) 



c=l c=l c,d=l 

n n 



If/) ki=Q = l-pY,9ca-2g aa + 2pg aa Y,9ca + (4 + l2h)g 2 aa + 



c=l c=l 

n n n 



+ ^ J29ca+Pj2 9cc9ca + P 2 J2 9dc9ca, (7) 



c=l c=l c,d=l 



where the notation 

1 d f d d k x d d k 2 



J 2 dk 2 J {k 2 + l){k 2 + + A: 2 + fc) 2 + i; 

]d u.jdl 



(8) 



fc 2 =0 



is used, and the redefinition g afe — > g ab / J is carried out. The diagram representation of 
the corresponding contributions to T^ 2 -*, rj^ and r^ 2,1 -* is shown in|2l 

However, the subsequent renormalization procedure for the vertex functions and the 
calculation of the (3 and 7 functions, which determine the renormalization group trans- 
formations for the interaction constants, are difficult due to the complicated structure 
of relations (JHJ) and (jlj) defining operations with matrices g ab . The step-like structure of 
the function g(x) established in |UE1IE] makes it possible to implement the renormaliza- 
tion procedure. In this paper, we restrict ourselves to the consideration of the one-step 
function q(x)\ 

9{x) = { (9) 
g u x < x < 1, 



a 



b 



4 /TV? 




2p 

#c6 ✓ & & • 



5W> 



-4 
6 b 
a a 




—4 fi'ab 



9ab 8 5W> < fi'ab I 9ab 



16 



6 6 
a 5U 6 



9ab > 9ab 



a b 

4 ~^a>^^j9bb 16 ~|^a 

a b CL dab b 



b b 

a 9ab b 



32 




9ab 



8 



8j9 



0ab 




— t — a 6> 

fi'aa I $ 9aa 8 

— ' a 6' 




- a 6> 

3aa y fi'aa 16 

- a 6- 
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0\ 0. 



2p 
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Figure 1: The diagram representation of contributions to the two-point (a), four- 
point T^' 1 ) (b), and with the inclusion (4>1) 2 (c) vertex functions in the one- and two-loop 
approximations with the corresponding weighting coefficients. 
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where the coordinate of the step < xo < 1 is an arbitrary parameter that does not 
evolve under scale transformations and remains the same as in the seed function go(x). 
As a result, the renormalization group transformations of the replica Hamiltonian with 
RSB are determined by the three parameters g,go,gi- 

The critical properties of the model can be revealed by analyzing the coefficients 
Px(tj,9o,9i) (* = 1>2,3), ^(g,g ,gi), and 7^(g,ft,9i) of the renormalization group 
Callan-Symanzik equation [Hj. We obtained the following expressions for the (3- and 
7 functions in the two-loop approximation in the form of series in the renormalized pa- 
rameters g,g and g^. 

Pi = -g + (8 + P )g 2 - P x g 2 ~ P (l-x ) gi 2 + [(8f-A0h + 20) P 

+ 16 / - 176 ft + 88]f + (24 ft - 8 / - 12) x P gg 2 + (24 ft - 8 / 

- 12) (1 - x ) P ~ggi - (16 ft - 8) x P g 3 - (16 ft - 8) (1 - x ) P gi\ 
P2 = -go + (4 + 2p) gg + (2 P x - 4) g 2 + 2 (1 - x )pg g 1 

+ [(8/-48/i + 28)p + 16/-48/i + 24]^ -[((32/i-16)a; (10) 

+ 8 - 32 h) P + 48 - 96 h]gg 2 - (32 h - 16) (1 - x )pgg gi 

+ [(48 h - 8 f - 20) x p - 32 h + 1% 3 + (32 ft - 8) (1 - x ) P g 2 gi 

+ (16/i-12-8/)(l-x )^i 2 , 

A = -(7i+px^ 2 -[p(x -2)+4]( 7l 2 + (4 + 2p)^ 1 + [(8/-48/i 

+ 28) P + 16 / - 48 h + 2A] 9l g 2 - (16 h - 8) x P gg 2 - [((8 - 16 h) x 

- 8)p + 48 - 96 h]u o9l 2 + (16 ft - 8) x ^ 3 + (8 ft - 8 / - 4) x P gig 2 
+ [(8 / - 24 h + 12) x p + (48 h - 8 / - 20) p + 16 - 32 ft] 

7^ = A(A-d)f(d)[( P + 2)g 2 - P x g 2 - P (l-x )g 1 2 ], 
1(f)2 = -( A - d )[( p + 2)g + p Xo g 0+p (l- Xo )g 1 -2{6h-2f-3){(p + 2)g 2 

- px go 2 -p (i - zo)^ 2 ) . 

In order to compare the results of this paper with those obtained in jH0 E], we reversed, by 
analogy with [H El IE] , the signs of the off-diagonal elements in the matrix: g a ^ — > —g a ^b- 
As a result, go and #i became positive definite. The integrals f(d) and h(d) were calculated 
numerically for 3 < d < 4. 

It is well known that the series used in perturbation theory are asymptotic, and vertices 
of the interaction of the order parameter fluctuations in the fluctuation domain r — > 



are sufficiently large to directly apply expressions (fTU|) . For this reason, in order to 
extract the physical information from those expressions, we use the generalized (for the 
three-parameter case) Pade-Borel method for summing asymptotic series. The direct and 
inverse Borel transformations have the form 

oo 

f(9,9o,9i) = £ Cij^g^gX = J e^F(gt, g t, git)dt, 

F(g,g ,gi) = £ T^f^g i glg\- 

In order to find the analytic continuation of the Borel image of a function, we use the 
following series in the auxiliary variable 9 

oo k k—i 

H~g, go, gi, #) = E^EE -^^-gWi ' '■ (12) 

k=0 i=0j=0 K - 

The Pade approximation [L/M] is applied to this series at the point 9 = 1. This technique 
was successfully used in || to describe the critical behavior of certain systems character- 
ized by several interaction vertices of the order parameter fluctuations. The symmetry 
conservation property of a system when applying a Pade approximant in the variable . 
becomes essential in the description of multivertex models. In this paper, we used the 
[2/1] approximant for the calculation of /^-functions in the two- loop approximation. 



3 Calculation results 

It is well known that the nature of the critical behavior is determined by the existence of 
a stable fixed point satisfying the system of equations 

Pi(g*,gZ,9l) = o (f = 1,2,3). (13) 

By numerically solving system ()13)) for functions found by the Pade Borel summation 
technique (for p = 1, 2, 3 ), three types of nontrivial fixed points were found in the physi- 
cally interesting domain of parameters g*,gQ,gl > (see Tables 1-3). The firsttype fixed 
point with g* ^ 0, g§ = g* = corresponds to the critical behavior of the homogeneous 
system; the second-type fixed point with g* ^ and g£ = g\ ^ corresponds to the 
critical behavior of the disordered system with replica symmetry; and the thirdtype fixed 
point with g* ^ 0, g$ = and g* ^ corresponds to the critical behavior of the disor- 
dered system with RSB. The values of and at the fixed point with RSB depend on the 
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coordinate of the step x . Tables 1-3 present the values of and for < x < 1 with the 
step Ax = 0, 1 The possibility to realize one or another type of the critical behavior for 
each p depends on the stability of the corresponding fixed point. The stability criterion 
of a fixed point reduces to a condition that the eigenvalues Aj of the matrix 

dPi{g*,g* Q ,gl) 

' h -i = % (14) 

belong to the complex right half-plane. 

An analysis of A, for every type of the fixed point (see Tables 1-3) provides the following 

conclusions. 

1) For the three-dimensional Ising model (p = 1), the second-type fixed point is stable 
(Table 1(a)). The complex values Ai and A2 for positive |Ai|, | A2 1 an d A3, show that 
the second-type fixed point, in contrast to the third-type one, is a stable focus in the 
parametric space (g,go,gi), and the renormalization group flows approach the second- 
type fixed point by a spiral curve. At the threshold dimension d c = 3.986 (see Table 1(b), 
(c)), the second-type fixed point looses stability (A3 changes sign). Since all other fixed 
points remain unstable in the entire range of the dimension variation 3 < d < 4, no critical 
behavior is realized in the system at 3.986 < d due to the replica symmetry breaking. 
The analysis of the behavior of renormalization group flows at 3.986 < d provides the 
following results. 

2) For the three-dimensional XY model (p — 2), small values of A« (see Table 2(a)) 
indicate that the second- type replica symmetric fixed point is weakly stable. However, 
already for the dimension d c — 3.1 (see Table 2(b), (c)), the third-type fixed point with 
RSB effects becomes stable. However, the critical behavior determined by this point is 
nonuniversal and depends on the parameter x and, therefore, on the concentration of 
impurities. A stability analysis of the third-type fixed point reveals that it is stable only 
in the interval < xq < x c (d), where threshold value of the parameter, which 
depends on the dimension of the system. For example, for d — 3.1 x c — 0.1; and, for 
d = 3.999 x c = 0.3 In the interval x c (d) < x < 1, all fixed points are unstable. 

3) For the isotropic three-dimensional Heisenberg model {p — 3), the first-type fixed 
point becomes stable (Table 3(a)), while at the other fixed points the constants g$ and gl 
take physically senseless negative values. Only at d c = 3.999, g$ and g\ take physically 
meaningful values for the third-type point, and this point becomes stable in the interval 
< x < 0.4 (Table 3(b), (c)). In the interval 0.4 < xq < 1, all fixed points are unstable. 



Table 1 

Coordinates of the FPs and eigenvalues of the stability matrix for p = 1: 

a) dimension d = 3.0 



type 


x 


9 


9o 


9i 


Ax A 2 


A3 


1 




U.l / f4 


U 


U 


U.OOoo — U.lbyz 


— u.ioyz 


2 




n 1 O A A 


0.0812 


0.0812 


0.5253 ± 0.0893? 


n oi 1 O 

U.zllz 


3 


0.0 


A 1 O A A 

0.1844 





0.0812 


r\ roro i n n o n o * 

0.5253 ± 0.0893? 


n n ono 

—0.0392 




0.1 


n i o a n 

0.1840 





r\ noon 

0.0829 


0.5352 ± 0.0983? 


n n /i no 

—0.0492 




r\ o 

0.2 


r\ i o o r - 

0.1835 





n n o a n 

0.0846 


0.5471 ± 0.1067? 


n n r nn 

—0.0599 




0.3 


n i o o a 

0.1830 





n n on o 

0.0863 


n tT/^nT i n 1 i oo ■ 

0.5607 ± 0.1133? 


n nTi 

—0.0712 




0.4 


n 1 oo /i 

0.1824 





n noon 

0.0880 


a r 'i r i nil on * 

0.5765 ± 0.1180? 


n n 

—0.0832 




0.5 


n 1 oi ^7 

0.1817 


n 




n n on r 

0.0895 


n rnn i n iono* 

0.5951 ± 0.1203? 


n nn r n 

—0.0959 




0.6 


0.1810 





r\ nn i n 

0.0910 


n 1 ^70 i n 1 1 on * 

0.6172 ± 0.1189? 


n 1 nn 

—0.1093 




0.7 


0.1802 





0.0924 


0.6439 ± 0.1114?- 


-0.1234 




0.8 


0.1793 





0.0936 


0.6760 ± 0.0921?- 


-0.1381 




0.9 


0.1784 





0.0947 


0.7135 ±0.0353?' 


-0.1534 




1.0 


0.1774 





0.0957 


0.8573 0.6536 


-0.1692 








b) dimension d 


= 3.985 




Type 


x 


9* 


9*o 


9*i 


Ai A2 


A3 


1 




0.0917 








0.6315 -0.4163 


-0.4163 


2 




0.1231 


0.1090 


0.1090 


0.6986 ±0.1311? 


0.0022 


3 


0.0 


0.1231 





0.1090 


0.7047 ±0.1069?' 


-0.0363 








c) dimension d 


= 3.986 




Type 


x 


9* 


9*o 


9*i 


Ai A2 


A3 


1 




0.0916 








0.6318 -0.4165 


-0.4165 


2 




0.1230 


0.1092 


0.1092 


0.6895 ±0.1453?' 


-0.0076 


3 


0.0 


0.1230 





0.1092 


0.7018 ±0.0935?' 


-0.0359 
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Table 2 

Coordinates of the FPs and eigenvalues of the stability matrix for p = 2: 

a) dimension d = 3.0 



Type 


X 


9* 


9o 


9*i 


Ai 


A 2 


A 3 


1 




0.155830 








0.667315 


-0.001672 


-0.001672 


2 




0.155831 


0.000584 


0.000584 


0.667312 


0.001682 


0.000004 


3 


0.0 


0.155831 





0.000584 


0.667313 


0.001683 


-0.000001 




0.1 


0.155831 





0.000614 


0.667313 


0.001684 


-0.000088 




0.2 


0.155831 





0.000648 


0.667313 


0.001685 


-0.000186 




0.3 


0.155831 





0.000686 


0.667313 


0.001686 


-0.000296 




0.4 


0.155831 





0.000729 


0.667313 


0.001687 


-0.000419 




0.5 


0.155831 





0.000778 


0.667313 


0.001687 


-0.000559 




0.6 


0.155831 





0.000833 


0.667313 


0.001688 


-0.000717 




0.7 


0.155831 





0.000896 


0.667314 


0.001690 


-0.000901 




0.8 


0.155831 





0.000971 


0.667314 


0.001692 


-0.001116 




0.9 


0.155831 





0.001058 


0.667315 


0.001694 


-0.001369 




1.0 


0.155830 





0.001163 


0.667316 


0.001696 


-0.001672 








b) dimension d 


= 3.10 






Type 




9* 


9*o 


9*i 


Ai 


A 2 


A 3 


1 




0.1499955 








0.689608 


-0.009539 


-0.009539 


2 




0.150017C 


0.00325 


0.00325 


0.689535 


0.009887 


-0.000003 


3 


0.0 


0.150017C 





0.00325 


0.689535 


0.009887 


0.000109 




0.1 


0.150016S 





0.00341 


0.689535 


0.009899 


-0.000401 




0.2 


0.1500167 





0.00360 


0.689536 


0.009926 


-0.000961 








c) dimension d - 


= 3.999 






Type 




9* 


9o 


9*i 


Ai 


A 2 


As 


1 




0.089762 








1.119442 


-0.133591 


-0.133591 


2 




0.092307 


0.036991 


0.036991 


1.103421 


0.227335 


-0.025378 


3 


0.0 


0.092307 





0.036991 


1.103421 


0.227335 


0.030783 




0.1 


0.092270 





0.038723 


1.102142 


0.235506 


0.021563 




0.2 


0.092205 





0.040559 


1.100913 


0.244667 


0.011135 




0.3 


0.092108 





0.042500 


1.099845 


0.254810 


-0.000648 




0.4 


0.091970 





0.044547 


1.099106 


0.265820 


-0.013939 
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Table 3 

Coordinates of the FPs and eigenvalues of the stability matrix for p = 3: 

a) dimension d = 3.0 



Type 

■J ± 


Xq 


9* 


9o 


9*i 


Ai 


A 2 


A 3 


1 




0.1383 








0.6814 


0.1315 


0.1315 


2 




0.1419 


-0.0359 


-0.0359 


0.6727 


-0.0891 


-0.1450 


3 


0.0 


0.1419 





-0.0359 


0.6727 


-0.0891 


-0.0058 




0.1 


0.1420 





-0.0382 


0.6727 


-0.0865 


0.0011 




0.2 


0.1420 





-0.0408 


0.6728 


-0.0836 


0.0088 




0.3 


0.1421 





-0.0439 


0.6730 


-0.0802 


0.0175 




0.4 


0.1420 





-0.0474 


0.6734 


-0.0764 


0.0273 




0.5 


0.1420 





-0.0516 


0.6738 


-0.0719 


0.0385 




0.6 


0.1418 





-0.0565 


0.6745 


-0.0668 


0.0515 




0.7 


0.1415 





-0.0625 


0.6755 


-0.0606 


0.0667 




0.8 


0.1409 





-0.0699 


0.6768 


-0.0533 


0.0845 




0.9 


0.1400 





-0.0793 


0.6787 


-0.0443 


0.1058 




1.0 


0.1383 





-0.0915 


0.6814 


-0.0331 


0.1315 



b) dimension d = 3.999 



Type 


x 


9* 


9o 


9*1 


Ai 


A 2 


A 3 


1 




0.0£ 


M91 


19 








1.113633 


-0.000820 


-0.000820 


2 




0.0£ 


$19* 


39 


0.000171 


0.000171 


1.113633 


0.000822 


-0.000228 


3 


0.0 


0.0£ 


$19* 


19 





0.000171 


1.113633 


0.000822 


0.000228 




0.1 


0.0£ 


$19* 


19 





0.000183 


1.113633 


0.000822 


0.000188 




0.2 


0.0£ 


$19* 


19 





0.000196 


1.113633 


0.000823 


0.000142 




0.3 


0.0£ 


$19* 


19 





0.000212 


1.113633 


0.000823 


0.000088 




0.4 


0.0£ 


\m 


19 





0.000230 


1.113633 


0.000823 


0.000025 




0.5 


0.0£ 


$19* 


19 





0.000251 


1.113633 


0.000823 


-0.000050 




0.6 


0.0£ 


\m 


19 





0.000277 


1.113633 


0.000824 


-0.000140 




0.7 


0.0£ 


$19* 


19 





0.000309 


1.113633 


0.000824 


-0.000251 




0.8 


0.0£ 


$19* 


19 





0.000350 


1.113633 


0.000825 


-0.000391 




0.9 


0.0£ 


\m 


19 





0.000402 


1.113633 


0.000826 


-0.000574 




1.0 


0.0£ 


\m 


19 





0.000473 


1.113633 


0.000828 


-0.000820 



12 



Table 4 

Critical exponents of the 3D models for RS FP's 



Model 


FP 




V 


V 


7 


P 


a 


Ising 


RSI 


this work 


0.028 


0.631 


1.244 


0.324 


0.107 






Ref. M 


0.031(4) 


0.630(2) 


1.241(2) 


0.325(2) 


0.110(5) 




RS2 


this work 


0.028 


0.672 


1.329 


0.345 


-0.015 






Ref. [HI 


0.030(3) 


0.678(10) 


1.330(17) 


0.349(5) 


-0.034(30) 


XY 


RSI 


this work 


0.029 


0.667 


1.318 


0.343 


-0.001 






Ref. M 


0.034(3) 


0.669(1) 


1.316(1) 


0.346(1) 


-0.007(6) 


Heisenberg 


RSI 


this work 


0.028 


0.697 


1.379 


0.369 


-0.092 






Ref. [IE! 


0.034(3) 


0.705(1) 


1.387(1) 


0.364(1) 


-0.115(9) 



Note that although the calculations indicate the stability of the impurity replica sym- 
metric second- type fixed point, there is reason to believe that, in the higher order approx- 
imations (as is the case for disordered systems considered without taking into account 
RSB effects [H]), the first-type fixed point, which corresponds to the critical behavior of 
the homogeneous system, will become stable. On the one hand, this is indicated by the 
very weak stability (A3 = 0.000004) of the second-type fixed point and by the fact that 
the threshold value of the order parameter p c = 2.0114 found in the two-loop approxi- 
mation, which separates the critical behavior domains determined by the first- (p > p c ) 
and second-type (p < p c ) fixed points, is very close to p — 2. This explains a very slow 
variation of the eigenvalues A; of the stability matrix for the disordered XY model with 
the variation of the system dimension (Table 2). On the other hand, the negative value 
of the critical heat capacity coefficient a of the XY model also suggests, according to 
the Harris criterion, that the critical behavior of the model is stable with respect to the 
influence of the quenched disorder and, therefore, it can be expected that p c < 2 in the 
higher order approximations. For example, the value p c = 1.912(4) was found in ^H] on 
the basis of the six-loop approximation with the use of the pseudo e-expansion and the 
Pade Borel Leroy summation technique with a thoroughly chosen fitting parameter. 

Because p c is very close to p = 2 for the XY model, one can expect that the calculations 
based on higher order approximations will substantially change the threshold dimension 
d c (p = 2), although for the Ising and Heisenberg models, the changes of d c {p) should 
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be small. This assumption is supported by the calculation of critical indexes for three- 
dimensional homogeneous models with p = 1,2,3 and the disordered Ising model. We 
performed these calculations in the two-loop approximation with the use of the Pade Borel 
summation technique (Table 4). The comparison of these results with the corresponding 
indexes reported in ^SJEh where the all-time accurate calculations for threedimensional 
models were performed in the six-loop approximation, shows that the difference in the 
values of critical indexes does not exceed 0.02. 

The values of the threshold dimensions d c (p), which separate the domain of critical 
behavior with RSB effects d c (p) < d < 4 from the domain where these effects are inessen- 
tial, can be considered as threshold dimensions that restrict the scope of the ^-expansion 
method as applied to the three-vertex model of the weakly disordered system and the 
corresponding results reported in ^EJ Ej. Our analysis also shows that the results of 
application of the e-expansion technique to multivertex statistical models are unreliable 
independently of the approximation order. This is explained by the competition between 
different types of fixed points in the parametric space of multivertex models, which usually 
does not allow one to pass to the limit as e — > 1 without crossing the marginal dimensions 
3 < d c < 4 separating the stability domains of different fixed points. 

To reveal the character of the behavior of a disordered system with RSB effects in the 
domain without stable critical states, we analyzed the phase portrait of the model based 
on the system of equations 

= Pi(g,go,gi), (15) 

which specifies phase trajectories in the space of vertices {g,go,gi)- An analysis shows 
(see Fig. 2) that, for the Ising model with d c = 3.986 at d > 3.986, where none of the 
fixed points is stable, the strong coupling regime with the renormalization group flows 
determined by (g, go, <?i)— 0^ 0) * s realized if g > g*. At the same time, at g < g*, the 
flows with (g, g , gi)— >(0, 0, 0) are realized, which asymptotically approach the Gauss fixed 
point (0,0,0) and then also tend to infinity along the axes g,go,gi- Such a behavior of 
the flows at g < g* is caused by the closeness of the system dimension d to four when the 
effect of fluctuations is negligibly small and the Gauss fixed point becomes an attractor. 
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Figure 2: The picture of renormalization group flows in the parametric space {g,go,9i) 
for the Ising model with the system dimension d = 3.99. 

4 Conclusions 

The renormalization group analysis of weakly disordered systems of an arbitrary dimen- 
sion in the range from three to four conducted in the two-loop approximation showed 
that the critical behavior of threedimensional systems is stable with respect to the effect 
of the replica symmetry breaking. In systems with a one-component order parameter, 
the critical behavior determined by the structural disorder with a replica symmetric fixed 
point is realized. The presence of weak disorder does not affect the critical behavior of 
multicomponent systems, although the proof of this fact for systems with p = 2 requires 
calculations with higher order approximations. 

Effects of the replica symmetry breaking manifest themselves only in disordered sys- 
tems with the dimension greater than three, and the threshold dimensions d c depend on 
the number of components of the order parameter p and the value of the parameter x - 
The predicted picture of the influence of replica symmetry breaking on the critical behav- 
ior of disordered systems with a dimension d > d c qualitatively agrees with the results 
reported in j3J El- The latter results were obtained on the basis of the e-expansion 
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technique. For systems with p = 1, RSB effects destroy the stable critical behavior, and 
the strong coupling regime is realized; for systems with p = 2 and 3, a domain of nonuni- 
versal critical behavior occurs at < xq < x c (d). For xq outside this interval, the system 
exhibits no critical behavior, as is the case at p = 1. 

The values of threshold dimensions d c {p): d c {p = 1) = 3.986, d c {p = 2) = 3.10, 
and d c (p = 3) = 3.999 which separate the domain of critical behavior with RSB effects 
d c (p) < d < 4) from the domain where these effects are insignificant, simultaneously 
specify the lower bound of the domain where the results obtained by e-expansion are 
applicable to the description of the model of weakly disordered systems with RSB effects 
jH 13 E] • R is noted that calculations carried out in higher order approximations of the 
theory can significantly change the threshold dimension dc for the XY model. On the 
other hand, changes in d c (p) for the Ising and Heisenberg models are expected to be 
small, which leaves the scope of the results obtained by the e-expansion technique close 
to dimension four. 

As the concentration of defects increases, one can expect a decrease of the threshold 
values d c down to d c < 3 beginning with a certain threshold concentration. In this case, 
the influence of replica symmetry breaking effects can be significant. Due to specific 
features of the manifestation of RSB effects, the concentration ns corresponding to the 
spin percolation threshold can play the role of the threshold concentration of defects for 
the Ising model, so that no stable critical behavior is observed at n > n s . For the XY and 
Heisenberg models, this role can be played by the concentration of defects corresponding 
to the impurity percolation threshold ni mp = 1 — n s with a nonuniversal critical behavior 
for n imp < n < n s and the absence of a stable critical behavior at n > n s . 
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